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ABSTRACT 

We have tested the two main theoretical models of bubbles around massive 
star clusters, Castor et al. and Chevalier & Clegg, against observations of the 
well studied Carina Nebula. The Castor et al. theory over-predicts the X-ray 
luminosity in the Carina bubble by a factor of 60 and expands too rapidly, by a 
factor of 4; if the correct radius and age are used, the predicted X-ray luminosity 
is even larger. In contrast, the Chevalier & Clegg model under-predicts the X- 
ray luminosity by a factor of 10. We modify the Castor et al. theory to take 
into account lower stellar wind mass loss rates, radiation pressure, gravity, and 
escape of or energy loss from the hot shocked gas. We argue that energy is 
advected rather than radiated from the bubble. We undertake a parameter study 
for reduced stellar mass loss rates and for various leakage rates and are able 
to find viable models. The X-ray surface brightness in Carina is highest close 
to the bubble wall, which is consistent with conductive evaporation from cold 
clouds. The picture that emerges is one in which the hot gas pressure is far 
below that found by dividing the time-integrated wind luminosity by the bubble 
volume; rather, the pressure in the hot gas is set by pressure equilibrium with 
the photoionized gas at T = 10'^ K. It follows that the shocked stellar winds are 
not dynamically important in forming the bubbles. 

Subject headings: ISM:bubbles 



INTRODUCTION 



Over a dyna mical time only 2% of the gas in a typical disk galaxy is turned into stars 
in spite of the fact that the disks are marginally gravitationally stable 
This puzzling result suggests that something other than gravity plays 
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a (negative) role in star formation. Passive support against collapse and subsequent star 
formation can come from magnetic fields or turbulence. Examination of star forming galaxies 
and individual star forming regions commonly shows evidence for more active support against 
gravity, for example, expanding bubbles. Active mechanisms, which come under the rubric 
of 'feedback', include energy or momentum input by active galactic nuclei, or from stars. 
Since most disk galaxies, including our own Milky Way, show little or no nuclear activity, 
and since bubbles in individual star forming regions are clearly not due to active galactic 
nuclei, in this paper we focus on feedback from stars. 

Stars of type B2 and earlier supply several types of feedback, including protostellar 
jets, main-sequence stellar winds, radiation pressure, gas pressure associated with ionizing 
radiation, and supernovae. The mechanical energy impar ted via a supernov a is comparable 
to that injected by the star's stellar wind over its lifetime (jCastor et al.lll975l ). However, it is 
clear that many observed bubbles have formed before any of the stars in the central cluster 
have exploded in a supernova. Since we are interested in bubble formation, we will ignore 
supernovae, leaving open the question of whether supernovae provide the feedback necessary 
to limit the rate of star formation. 

There are two co mpeting theoretical mod els of stellar wind b ubbles in th e literature: 



that of Castor et al. (i.e. Castor et al. 



19751) andlWeaver et all (119771 )) and that of lChevaher fc Cleee 



(1l985l ). as implemented in, e.g., IStevens &: Hartwelll (l2003l ). In the theory explored by Cas- 
tor, the stellar-wind shocked gas is confined by a cool shell of swept up interstellar medium 
(ISM). In contrast, the Chevalier & Clegg theory ignores any surrounding material and sim- 
ply has a steady-state wind. The temperatures, pressures and predicted sizes of the X-ray 
emitting regions in the two theories are significantly different. 

However, these differences can be masked if the stellar content of the bubble is not well 



constrained. Recent work by ISmithI (l2006al ) provides a detailed accounting of the number 
of early type stars in the well studied Carina nebula. Combined with the wealth of multi- 
wavelength observations available for Carina, this allows for a stringent test of the models. 

Observations of star forming regions in the Milky Way show that the gas density is highly 
variable. In particular, the projected surface den si ty is consist e nt wit h a log-normal distri- 
bution, e.g., ICoodman. Pineda, fc Schned (120081 ): IWong et al.l (120081 ). These observations 
are consistent with both analytic and numerical studies of supe rsonic turbulent flows, which 
predict log-normal density and colurn n density distributions (jPassot fc Vazquez-Semadeni 
19981 : lOstriker Stone fc Gammid l200ll ). These results call into question one of the funda- 
mental assumptions of the Castor et al. model, that the gas surrounding the central star 
or star cluster has a rather uniform density distribution. If the mass distribution is non- 
uniform, the shell swept up by the shocked stellar wind is likely to have holes. If so, this will 
lead to incomplete conflnement of the hot gas and a consequent reduction in the pressure and 



- 3 - 



associated X-ray emission of the hot gas. Such a situation would be intermediate between 
the Castor et ah model of a completely confined wind and the Chevalier & Clegg picture of 
a free-flowing wind. We explore such intermediate models below. 

This paper is organized as follows. In the next section we summarize some relevant 
observations of Carina, as well as some order of magnitude estimates of the pressures and 
forces acting on the surrounding cold, dense gas. In ^ (and the appendices) we describe 
the Chevalier & Clegg model, the Castor et al. model, and our modifications to the Castor 
et al. models. We also describe our numerical methods. In §l]we present the results of our 
modehng, including a parameter study of the (modified) Castor et al. bubbles. In ^ we 
discuss the X-ray emission. In §H] we discuss the results and compare to previous work. We 
give our conclusions in the final section. 



2. OBSERVATIONS OF CARINA 



The distance to Cari na is 2.3 kpc (lAUen fc Hillierl Il993l : ISmithI l2006bl ). The complex 
contains about 70 O stars (ISmithll2006al ). of which 47 (including evolved massive stars) reside 
in Trumpler 16 (Tr 16), a star cluster with a radius r ~ 3pc. This makes Carina one of the 
most massive (and best studied) star formation regions in the Milky Way. The bolometric 
luminosity emerging from t he stars is 9.6 x 10^°erg s~^. This is about twice the far infrared 
luminosity of the complex (jSmith fc Brooks! l2007al ). indicating that half the stellar flux is 
intercepted by dust grains within ~ 10 — 20 pc of the stars. 

We note that this rather low ratio of infrared to bolometric luminosity is consistent with 
the fact that we can see Tr 16 in the optical, and with turbulent models of the ISM, but it 
is inconsistent with the assumption of the Castor et al. model that the ISM and swept-up 
bubble shell surrounding the cluster is homogeneous. It suggests that hot gas can escape the 
bubble, a point we return to later. 

The number of ionizing photons emitted per second is Q = 9 x 10^'^ s^^. The kinetic 
energy of the stellar winds is ~ 3.5 x 10'^^ erg s~^, assuming standard estimates of the 
mass loss rates; we argue below that this is likely to be an overestimate by a factor of 3 — 10. 
The massive stars of Tr 16 are responsible for about 70% of the bolometric and stellar wind 
luminosity and of the ionizing flux of the entire nebula; see Table [TJ 

A reasonably reliable estimate for the age of Tr 16 is given by the age of rj Carina, 
currently the most massive star in the region, t] Carina is ~ 3.6 million years old, since it has 
clearly evolved off the main seq uence, while estimat es of its initial mass based on its current 
luminosity and stellar models (IBressan et al.l Il993l ) are in the range of 12OM0. We thus 
assume the cluster is approximately 3.6 million years old. It is believed that no supernovae 



-4- 



have occurred in Tr 16 flRetallacklll983l : iLopez fc Meaburnlll984l : IWhiteoaklll994l ). but in any 

case the integrated energy deposited by the winds, 4 x 10^^ erg, substantially exceeds the 
energy input of a supernova, of order 10^^ erg. 



For the 17th anniversary of the Hubble Space Telescope, ISmith fc Brooks! (l2007bl ) took 
a 50 light-year-wide mosaic of the centre of the Carina nebula. This image has unprecedented 
detail of globules within the bubble and star birth at the edge of the bubble. Triggered star 
formation around the edge of the bubble in Carina is also observed in X-ray and Near-IR 
studies of Carina fISanchawala et all 120071 ). There are many directional markers, such as 
pillars, showing that the bubble edge is dominated by the effects of Tr 16. The images show 
that the main bubble shel l is markedly non-spheri cal. Smaller bubble structures within the 
main bubble can be seen. ISmith fc Brookj (j2007al ) compiled a multi-wavelength analysis of 
the Carina nebula enabling approximate edges to the bubble to be located. We estimate 
radii in the range 10 — 20 pc. 

There have been many observatio ns of Carina in X-rays, which establish the existence of 



diffuse emission from the nebula itself ( Seward et al.lll979l : lFlaccomidl2005l: Hamaguchi et al. 



2007al lb[). The largest area coverage of the Carina region is provided by ISeward et al.l (119791 ) , 
who used the Einstein Observatory to find a diffuse X-ray emission of 4.8 x lO^'^erg s~^ for 
the 0.5 — 3.0 keV band (scaled from their distance estimate of 2.6 kpc to 2.3 kpc). 

The pressure inside the bubble surrounding Tr 16 contributes directly to the outward 
force on the surrounding material. We estimate the pressure in several ways. 

First, the pressur e can be calculated from the free- free radio emission, e.g., as reported 



m 



Huchtmeier fc Day! (1l975h . who measured a flux of Si, = 1.43 x 10^ Jy. From calculations 



of the expected Bremsstrahlung radio flux at the observer 



uhii = 3.6 — 
\pc 



-) 
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e^i'i cm 



(1) 



where enii is the volume filling factor of the H II gas inside the bubble. For the radius of 
the bubble, i? = lOpc and distance, D = 2.3 kpc, we find uhii = 100 cm~^ for Thjj = lO'^K. 
Thus, P = 1.4 X 10-i°e^y/ dynes cm'^. 

As a consistency check, we note that pressure of an H II region powered by a cluster 
emitting Q ionizing photons per second is given by 



HII 



3Q 



Airr^areHii 



hTnii 



where ~ 



4 X 10 ^^cm^s ^ is the recombination coefficient, and T^n 
Q = 9 X 10^° s~^ and r ^ 10 pc, we find Phii ~ 2 x 10~^°e^y^ dynes cm~^. 



(2) 



10^ K. Using 
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Second, we can use the ob served X-ray flux to estimate the pressure in the hot (Tj 



10^ K) gas. From Figure 1 in ISeward et ahl (Il979l ) we estimate the radius of the diffuse 



X-ray emission region to be ~ 12 pc, similar to the size of the free-free emission region. 
Combined with the diffuse X-ray luminosity and T^. = 6 x 10^ K, this yields n^^^ ^ 9 x 10~^ 
and ~ 7 X 10~^^ dynes cm~^, which we take to be equal to the pressure of the radio 
emitting gas. 

The surface brightness of both the radio and X-ray emission is not uniform, indicating 
that the pressure is not constant inside the bubble. For example, for the two classical and 
smaller scale embedded H II regions Carina I and Carina II, using the fluxes and sizes from 



Huchtmeier fc Day! (Il975l ) we flnd a pressure 9 x 10 dynes cm . Similarly, ISeward et al. 



(1l979h estimate the density of the hot X-ray emitting gas in the prominent r ~ 1.25 pc diffuse 



patches to be n ~ 0.4, giving P ^ 3 x 10 dynes cm ^. 

Finally, we can estimate the pressure from the outflow seen to be emerging from the 
globule known as 'the flnger', which has a dense ionization front and photoevaporative flow, 
with a terminal shock resolved close to the edge. Estimates of the ionization front density 
can be obtained from the S II emission line, and from the Ha emission measure. These 
give electron densities of 2000-6000 cm~'^ at various positions. The dense ionization front 
drives a photoevaporative flow out into the H II region until it reaches a terminal shock 
where the pressure is balanced with the ambient medium. To calculate Pjp assume the 
temperature is 10^ K due to hydrogen ionization and the flow is spherical. Nathan Smith 
(private communication) made two estimates, one using the tip of the flnger with a very 
small ionization front radius but a relatively high number density of 6000 cm~^ and a second 
using the main part of the flst with a larger ionization front radius and a lower density of 
3300 cm~^. These yield 5.4 x 10^^'^dyne cm~^ and 7.3 x 10~^'^dyne cm~^, respectively. 

These rather direct measurements of the pressure in several gas phases can be compared 
to two other pressures, the pressure exerted by the self-gravity of the surrounding molecular 
gas (as estimated from CO measurements) and the pressure of radiatiori from the stars 



in Tr 16 on dust grains in the interstellar medium. lYonekura et al.l (120051 ) flnd a mass of 



molecular gas M„ ^ 3.5 x IO^Mq in a region of radius r ^ 20 pc, yielding a surface density 



'9 

E « 6 X 10~^gcm~^. Assuming this gas is in hydrostatic equilibrium (there is evidence 
that it is expanding outward, but we ignore this for the moment) the dynamical pressure is 
P ^ TiGT? ^8x 10-1° dynes cm-2. 

The near equality between the pressure in the X-ray gas and that of the H II gas likely 
reflects the fact that the sound crossing time (of the cooler gas) is similar to the dynamical 
time or the age of the system. Like the hot gas, the H II gas will 'head for the exits'; pressure 
gradients will drive the warm gas toward any holes in the bubble shell, reducing the pressure 
on a (H II gas) sound crossing time. 
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Unlike the hot gas, the velocity of the warm gas is comparable to the velocity of the 
bubble wall, so a pressure driven flow will not reduce the H II pressure in the bubble by a 
large (factor of 10 or more) amount. Thus it appears natural that the H II gas will have a 
pressure similar to that estimated from the ionizing photon luminosity. Conversely, the hot 
gas pressure will, if it ever reaches values similar to L^r/V , drive outflows through holes in 
the bubble that have velocities of order the sound speed of the hot gas, which is two orders 
of magnitude larger than the velocity of the bubble wall. This will lower the pressure of the 
hot gas until it approaches that of the H II gas, at which point the pressure gradients in the 
two components are similar. 

Using the dynamical argument that the H II pressure and the X-ray gas pressure should 
be roughly equal, we can draw one more conclusion: that enii is of order unity. To do so, 
set the right hand side of eqn. [2] equal to the (observed) pressure of the x-ray gas and solve 
for euii, then use the known values of Q, r, and Thii- 

The pressure exerted by radiation is Prad = -^bo«/(4vrr^c) ~ 4xl0~^''(10pc/r)^ dynes cm~^. 
This near-equality, unlike that between the H II and X-ray gas pressure, is a bit of a coinci- 
dence. The radiation pressure scales as while the H II gas pressure scales as i?"^/^; the 
radiation pressure was more important in the past, while the gas pressure will become more 
important as the bubble expands in the future. 

All of these observationally driven estimates are well below the pressure estimated by 
dividing the total stellar wind kinetic energy by the volume of the bubble; the latter is 
P-wind = 3.6 X 10~''(10 pc/r)^(Lu;/4 x 10^* erg s~^) dynes cm~^. Either the wind has expanded 
over a volume larger by a factor of ~ 1000, the wind luminosity is overestimated by a similar 
factor, the shocked gas loses most of its internal energy (but not by radiation in the X-ray 
band), or some combination of all three. We note that if the hot gas expanded at its sound 
speed, it would reach radii of 1 kpc, so the notion that the hot gas has expanded through 
holes in the ISM to occupy region with a typical size of 100 pc is not implausible. 

Assuming an age of 3.6 Myrs, we find an average expansion veloci ty for the bubble wall o f 
~ ±3 (-R/lOpc) km s~^, similar to the 5 km s~^ splitting seen in CO ( Grabelsky et al.lll988l ). 
Observations of ionized gas typically find larger velocity splittings, of order ±20 km s~ \ 
ranging up to several hundred kilometers per second e.g. IWalborn fc Hesserl (1 19751 . Il982l ). 
but we regard this as a measure of the velocity of gas in the interior of the bubble, not a 
measurement of the velocity of the bubble wall. 

To sum up, the warm and hot gas appear to be in pressure equilibrium, with a pressure 
Pg ^ 2 X 10~^° dynes cm~^. This is similar to the radiation pressure, and to the pressure 
exerted by the surrounding, self-gravitating molecular cloud. The kinetic energy input from 
stellar winds greatly exceeds the integral of this pressure over the apparent size of the H II 
region, R ^ 10 pc, suggesting that the winds have escaped to occupy a much larger volume. 
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or been otherwise dissipated. 



3. METHODS 



3.1. Chevalier and Clegg bubble 



Stevens fc Hartwelll (120031 ) applied the IChevalier &: Clegd (Il985l ) steady state wind to 



super star clusters. The flow is assumed to pass through the sonic point at the cluster 
edge where the source terms become zero. The equations describing the flow are given in 
Appendix I. 

Numerical integration of the flow (equations (1A1|) to (]A3|) ) using iterative calculations 
of the mach number (equation (]A6|) or (]A7I) ) at every step enables the flow to be accu- 
rately calculated. The X-ray luminosity from Bremsstrahlung is numerically integrated 
from the resulting temperature and density proflles. The actual X-ray luminosity will be 



somewhat larger than this estimate by about a factor of three at 
Sutherland fc Dopital Jl993h . 



10'' K, see, e.g., 



3.2. Castor bubble 

In the Castor et al. theory the stars are assumed to blow a constant, spherically sym- 
metric stellar wind that interacts with the ambient ISM and produces a bubble. For the 
flrst few hundred years the stellar wind freely expands at the wind velocity. Once the wind 
shocks, there is a period of adiabatic expansion for a few thousand years. The bubble then 
enters the "snow plough" stage when swept up mass densities reach the critical point for 
radiative cooling and a thin shell forms. In the Castor scenario the snow plough stage is the 
most readily observed phase as it is by far the longest lived. During this time the bubble 
has a four zone structure, (Fig. [H^i)): a) A hypersonic stellar wind, b) A hot (T ~ 10^ K), 
almost isobaric region consisting of shocked stellar wind mixed with a small fraction of the 
swept-up interstellar gas. c) A thin, dense, cold (T << 10^ K) shell containing most of the 
swept up material at Rb. d) Ambient interstellar gas. 

One of the most important features of the snow plough structure is the coohng between 
zones b and c. The shell, c, is dense and cold due to efficient radiative coohng. However, 
zone b is hot as its density is too low to radiatively cool. Thermal conduction between b and 
c causes gas from the inner edge of c to evaporate into b. Although the mass loss from c is 
negligible compared to the accretion of mass from swept up ISM, it dominates over stellar 
wind as a source of matter in b. The bubble shell is pushed outward by the pressure from 
the confined hot gas. The equations describing the evolution are given in Appendix II. 
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Unlike the Chevalier & Clegg theory, the Castor theory takes into consideration the mass 
swept up frora the surrounding nebula, r e quiring an ISM or GMC m odel. The original papers 
( ICastor et al.l 1 19751 : IWeaver et al.l 119771 : IChevalier &: Clegd Il985l ) assumed a homogeneous 
ISM. However, in our models for Carina we sometimes employ an isothermal sphere as an 
approximation for the GMC surrounding the cluster. The values for t he cluster gas r nass and 
radius, and the surrounding GMC mass and radius were taken from [Murray et al.l (jin prepi ) 
assuming a density profile in the cluster of and in the GMC of r~^. The 'standard' ISM 
model used for our bubble models was Rd = 1 pc, = 2.6 x IO^Mq, Rq = 24.62 pc, and 
Mg = 3.2 X IO^Mq, where M^i is the mass in gas and dust remaining in the cluster after 
star formation assuming a star formation efficiency of ~ 35%. The Rq was chosen so that 
the density at the boundary between the cluster and GMC was continuous. 

In other cases we assume that the density falls off as p(r) ~ 1/r, or even that the density 
is (on average) uniform. 

In the Castor models we assume that the cluster produces a point source wind from the 
centre of the cluster. Although t his does not accurately represent the stellar wind within 
the cluster radius, simulations by ICanto et al.l (120001 ) show that the approximation is valid 
once outside the cluster radius. 

To model the Castor bubble, we used the integrate_ode program from Numerical Recipes 
(jPress et al.l Il996l ) to simultaneously integrate the momentum, radius, mas s and energy 



evolut ion equations (equations ( IBllI) to ( 1B14I) ) from the initial conditions in IWeaver et al. 
Jl977h (equations flBTSj) to flB23|) l 



The X-ray luminosity was calcula ted taking into acc ount the radial profile of density 
and temperature through the bubble (IWeaver et al.l 119771 ). When the X-ray luminosity is 
calculated the bubble is split into 5000 shells of equal width between Ra and Rb and the 
X-ray luminosity in the Einstein band from Bremsstrahlung radiation is calculated for each 
shell. A resolution check was conducted at many different bubble radii to ensure accuracy. 



3.2.1. Modifications to the Castor bubble model - stellar mass loss rates 



The generally accepted mass loss rate of O stars assumes a smooth (not clumpy) out- 
flow throughout the main sequence lifetime. However, recent observations show that stellar 
winds are clumpy and hence the rnass loss rates are likely lower than previ ously accepted 
JPuUerton et allbooi ISouret et allbood : lEvans et al. l 120041 : IPuls et allEooel l These obser- 



vations suggests that the mass loss rate is at least three times and possibly more than 30 
times lower than previously expected. We introduce a to parametrize the reduction in the 
the mass loss rate: = aM^,, original and = a L-^ ^original- 
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3.2.2. Radiation pressure 



The far infrared luminosity in the direction of Carina is roughly half the total stellar 
luminosity. It follows that half or more of the momentum carried by light emitted by the 
stars is deposited into the gas near the stars. Thus, the radiation pressure contributes to 
the outward force on the bubble shell (where the UV optical depth is large) 

- Lbol 



rad 



V- 



(3) 



where rj is of order 1/2, c is the speed of light and Lboi is the total bolometric luminosity 
from the cluster. 



Using a mult i- wavelength analysis of the Carina Nebula ISmith fc Brooks! (l2007ar) es- 
timate dust temperatures and calculate t he total mass in dust. Smith fc Brooksl ~ 2007a ) 
calculat e the masses using t he models of iGilmanI (119741 ). We employed more recent dust 
models (iLaor fc Draindll993l ): the mass in dust can be calculated using the mass per grain, 
47ra^p/3 and luminosity per grain = Q^aT^Aira^ 

ap 



(4) 



where a = effective grain radius, a is the Stefan-Boltzmann constant, = mean thermal 
emissivity, = dust temperature, and = the dust luminosity. No silicate features are seen 
so we can assume the grains are carbon. We also assume that gra ins are small (a < 0.2u m). 
The Qe values for each temperature can be read from the graph in lLaor fc Draind ( 19931 ). In 



0) 



particular, Qp ^ 3 x 10~ '^ for T^^ = 35 K. This is factor of 9 larger than the value employed by 
Smith &: Brooks! (l2007a! ): conse quently, we find a total du st mass ~ 1000(p/3gcm~^)M( 
Table [21 a factor 9 smaller than ISmith fc Brooks! (l2007a! ) . 

Assuming a dust to gas ratio of 1:100 the total mass for t he surrounding n ebula is 



IO^Mq, corresponding to an Ay ~ 1 (20 pc/i?;,)^. In contrast, ISmith &: Brooka (l2007a! ) 



find 9.6 x lO^M©, corresponding to an Ay ~ 10 (20pc/-R6 



Smith fc Brooks! (l2007a! ) find To ^ 35 K, which we argue supports the lower mass 



estimate, as follows. Using the expression for dust temperature 



Td = 49/ 



-l/(4+/3) 



2 X 10 



17 



cm~^ 



2/(4+/3) 



^bol 



\ IO^Lq 



l/{4+/3) 



(5) 



from IScoville &: Kwan! (1l976l ) with / = Qabs{^ = 50fi) ~ 0.01 we find that the typical 
distance between a dust grain and Tr 16 is ~ 24 pc (taking (3 = 1). We note that this 
is similar to the estimated bubble size. Since the bulk of the stellar radiation is initially 
emitted in the UV, it will be absorbed at < 1, consistent with the lower mass estimate 
given above. 
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3.2.3. Self- Gravity 

As the shell sweeps up a substantial amount of mass, the self-gravity of the shell and 
gravity between the shell and the stars and interior gas should be considered in the momen- 
tum equation. 

_ _GMl _ GM,{M, + M,t,rs) ... 



3.2.4. Leakage of hot gas 

The Castor theory assumes expansion through an homogeneous ISM. However, as noted 
in the introduction, turbulence is expected to produce density fluctuations in the ISM, while 
observations show that the column density is log-normally distributed. Thus, sections of the 
bubble shell with low column will expand faster and high column sections will expand more 
slowly. This uneven expansion will result in gaps in the shell, allowing the hot gas to escape 
from the bubble. 

We investigated the effects of holes in the expanding shell using a toy model. Our model 
consisted of a Castor-style bubble with a shell cover fraction Cfi when Cf = 1 the shell has 
no holes and when Cf = the shell is covered with holes, i.e. no shell exists. These holes 
allow the interior gas to escape at the speed of sound, Cg = {5ki,T/{3mp)y^'^, with a mass 
flux and energy flux given by 

M, = -(l-C/)47ri?,V., (7) 

E, = -il-Cf)A7TRl^p,cl (8) 

where Rb is the radius of the bubble interior, pb is the density inside the shell, which is 
assumed to be homogeneous and 7 = 5/3 throught this paper. 

We are not concerned with the bubble evolution within the cluster radius (of order 3 pc 
for Tr 16). Thus, we neglect gravity and radiation pressure for r < Rd] once the bubble 
expands, under the influence of the other forces, beyond this radius, we smoothly turn on 
gravity and radiation pressure. 

With hot gas from the interior able to leak away, it is important to consider the location 
of the stellar wind shock, Ra (see Fig. [1]): if -Ra approaches the shell radius at any time our 
model breaks down. We use the assumed ambient ISM density profile for the regions where 
the shell is still intact and negligible density where there are holes. Although this is not an 
accurate description of a clumpy ISM, it does allow the shell to leak while still gaining mass 
as it expands through a variable density ISM. The gas and radiation pressure only acts on 
the remaining shell (i.e. area = CfATcRf ). 
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Without leaking the radiative losses in region b were negligible. However, with leaking 
this may not always be the case. To ensure these losses are taken into account when sig- 
nificant we calculate the radiative luminosity of region b, Lf,, and include it in the energy 
equation. 

The sizes of the holes in the shell may vary as the bubble expands, so that Cf = Cf{t) 
or Cf{r). To investigate the effects of different hole expansion models, we conducted runs 
with many different power law relationships with time (including Cf{t) = const.), which we 
discuss below. 



3.2.5. Final bubble model 

Considering radiation pressure, gas pressure, stellar winds, gravity and leaking the evo- 
lution equations for momentum, radius, mass and energy become 

— = 47ri?,n + — , (9) 



dRb Vc 



dt Mc' 

dMb _ 5/2d2/d d \-i ^n~^l^ t ^^ \a^d2 



(10) 



C,rj'Ri{Rb - RaV - ^^^Lb - (1 - Cf)ATiRipbCs, (11) 
^ = L^-A^RlPb^-^-Lb-\{l-CMr.Rlpbcl. (12) 

respectively, see Fig. [l](n). 

Vc is the momentum of the shell, t is time, Rb is the radius of the bubble, Pb is the gas 
pressure in the bubble interior, Lboi is the total bolometric luminosity of the stars, c is the 
speed of light, G is the gravitational constant, M^. is the mass of gas in the bubble shell, Mb 
is the mass of gas in the bubble interior, = lO^M^u^ is the mass of stars in the cluster, 
Ci is a combination of constants associated with conduction, Tb is the temperature of the 
gas in the bubble interior, Ra is the radius of the stellar wind shock, kb is the Boltzmann 
constant. Lb is the luminosity emitted from the bubble interior, is the covering fraction of 
the shell, pb is the mass density of the bubble interior, is the speed of sound of the bubble 
interior, Eb is the energy in the gas in the bubble interior, is the stellar wind luminosity. 

The terms on the right hand side of the momentum equation ([9]) represent, in order, the 
outward force from the gas pressure, the outward force from radiation, the inward force from 
the self gravity of the bubble shell, and the inward force from gravity between the star cluster 
(with a minor contribution from the gas in the bubble interior) and the bubble shell. The 
terms on the right hand side of the mass equation f|TT]) are the mass input associated with 
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conduction evaporating the bubble wall, a (downward) correction in the conductive mass 
loss rate associated with radiative cooling, and mass loss from hot gas escaping the bubble. 
The terms on the right hand side of the energy equation ([T^ are stellar wind luminosity, 
adiabatic expansion, radiative cooling, and energy advection associated with the loss of hot 
gas through the bubble wall. 



4. RESULTS 

4.1. Original stellar wind models 

The Chevalier & Clegg model predicts a free-free X-ray luminosity in the Einstein band 
(0.5 — 3keV) of L.j. = 3.75 x 10^^ erg s~^ for Carina. The Castor model predicts an X-ray 
luminosity of = 5.4 x 10'^^ erg s"^ for Rb = 13 pc and = 3.8 x 10^^ erg s~^ for Rb = 20 pc. 
Plots of the models can be seen in Fig. [2] and Fig. [31 Both of these models considered X-ray 
emission from r > R^ to match the observations, in which the cluster region is ignored due 
to contaminating stellar and wind-wind coUisional X-ray emission. 

In the Chevalier & Clegg model the inner region of the wind dominates the X-ray 
luminosity so a strong radial gradient in the X-ray surface brightness is predicted; beyond 
a few parsecs the X-ray luminosity is negligible. In contrast, in the Castor et al. model the 
density is almost homogeneous within the bubble, so very little radial surface brightness 
gradient is expected, aside from a sharp edge at Rb- 



Observations of diffuse X-ray emission in Carina give = 6. 1 x 10^^ erg s ^ (ISeward et al 



19791 ): those authors assumed the distance to the cluster was 2.6 kpc. If we assume the dis- 
tance is 2.3 kpc, the luminosity becomes 4.8 x lO^^'erg s"^. The observed X-ray emission is 
ten times higher than predicted by Chevalier & Clegg model but one hundred times lower 
than predicted by the Castor model. Fig. HI The observed pressure, as estimated by 'the 
finger' and radio emission, is closer to the calculated Castor pressure when the bubble has 
the currently observed radius Rb ~ 20 pc. 

As noted elsewhere, the Chevalier & Clegg model predicts a pressure that has a strong 
radial gradient, so the pressure is dependent upon the location within the bubble rather than 
the bubble size. However, the Chevalier & Clegg pressure is far below the observed pressure 
at most points in the bubble. 

It is unclear how long it takes the bubble to reach the cluster radius, as the expansion 
time will depend upon the rather complicated wind-wind interactions within the cluster. To 
properly understand the early evolution of bubbles, full hydrodynamical models are needed. 
However, for now we will assume it is much less than a million years, as material will be 
displaced during star formation due to pre-main sequence jets and outflows. We assume 
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that the majority of the time taken to expand to the observed size is accumulated when 
Rci < r < Robs- Using this assumption we found that to expand to 20 parsecs the Castor 
model took only 0.90 million years, far shorter than the estimated age of Carina (around 3.6 
Myrs). 

It follows that the situation for the Castor et al. model is actually worse than the 
impression conveyed by the Figure, for the following reason. The actual age of the bubble is a 
factor of four higher than in the model, so the energy content of the bubble is underestimated 
in the model by the same factor, as is the pressure. Given the weak dependence of the 
temperature on time, the X-ray luminosity scales as the square of the pressure; if the radius, 
wind luminosity, and age of the bubble are set by observations, then the X-ray luminosity 
would be a factor of 16 larger than in the model. 



4.2. Modified Castor et al. model; leakage 

Taken together, the overestimate of the pressure and X-ray flux, and the underestimate 
of the bubble age (or overestimate of the bubble size at a fixed age) suggest that the outward 
force in the Castor et al. model is overestimated. 

The inclusion of radiation pressure and gravity for the Tr 16 model makes only a small 
difference to the evolution of the bubble. For Tr 16 \FHn\ ~ \Frp\ ~ |-Fgrai)|; in the model, 
unlike the case in Carina, all are much smaller than the gas-pressure force exerted by the 
shocked stellar wind, assuming no leakage and a standard O star mass-loss rate. 

In an attempt to produce models in closer agreement with observations, it is clear that 
the internal pressure of the bubble must be reduced. One way to do this is to allow hot gas 
to leak out of the bubble (C/ < 1), or equivalently to allow energy to leak out of the hot gas, 
via radiation or conduction; another is to reduce the mass loss rate and hence the kinetic 
luminosity of the stellar winds {a < 1). 

Observation of massive clusters support the possibility of leaking. Firstly, for Carina 
and many other cluster s you can see the stars, s howing that there must be a gap in the shell 



along our li ne of sight (jSniith fc Brooks! l2007bl ). Secondly, recent observations of Ml 7 and 
RCW49 by Povich et al. ( 20081 ) show stellar wind bow shocks around O stars at the edge 



and outside of the bubbles suggestive of large scale gas outflow from the H II region (see 
their Figure 2, especially RCW49-S1). 

Our initial simulations of a leaking bubble looked at three examples in detail: Cf = 1.0, 
Cf = 0.65, and Cf = 0.3. The results are given in Fig. O Not surprisingly, these simulations 
reveal that as C/ decreases the shell expands more slowly, the bubble interior pressure is 
lower, as is the predicted X-ray luminosity. 
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The effects of varying Cf witli time were tlien investigated. Tliese simulations started 
witfi Cf = 0.6 at t = 10'^ years and evolved Cf{t) witfi different power laws until t = 3Myrs, 
when Cf = 0.2. The simulations show that expansion rate, pressure, temperature and 
density differ very little with different power-law exponents; however, the X-ray luminosity 
does depend on the exponent. The difference in X-ray luminosity is due to the combination 
of the small changes in the parameters giving larger X-ray luminosities for larger exponents; 
larger exponents imply Cf being large for longer times. Thus, less material escapes over the 
lifetime of the bubble: Mescape oc {1 — Cf) and the X-ray luminosity is higher during the 
expansion when the power-law exponent is larger. 

To see the way Cf and the run of density in the ISM affected the evolution of the 
bubble we ran four different models. We considered Cf = 1.0 or 0.5 and either isothermal or 
constant density models (with the same total mass). Fig. [61 From this study we find that a 
homogeneous medium leads to a faster expansion; the mass swept up in the bubble shell is 
less at a given radius in the constant density model. However, the isothermal sphere model 
does approach the homogeneous model near the end of the simulations when both reach Rq 
and hence the same shell mass. The models with Cf = 0.5 expanded more slowly and had 
significantly lower X-ray luminosities than models with Cf = 1.0, as expected. 

We conducted a parameter study in a and C/ to see if any combination would result in 
acceptable fits to the observed bubble in Carina. Models were calculated for 0.15 < Cf < 0.99 
and 0.01 < a < 1.00, stepping both parameters at intervals of 0.01. The results are shown in 
Fig. [71 contours are only plotted for Ri, < 50 pc at a bubble age of 3 Myrs. This parameter 
space study was run for four different ISM models. Runs shown in part (a) of the Figure are 
for an isothermal sphere model with Rq = 24 pc; run (b) shows an constant density model 
with Rg = 24 pc; run (c) is for an isothermal sphere model with Rq = 60 pc; and run (d) is 
for a constant density model with Rq = 60 pc. 

The results show that isothermal sphere models gives smaller radii and slightly higher 
X-ray luminosity for any given {a,Cf) pair. The final radius at 3 million years increases 
with increasing a and with increasing Cf, as expected. For the isothermal ISM model the 
shell actually collapses inwards for low values of a and low values of Cf, a result of the self 
gravity of the shell exerting an inward force larger than the sum of the gas pressure force 
and radiation pressure force on the shell. 

The X-ray luminosity increases with increasing a and increasing Cf, again as expected. 
The region of parameter space that best fits the observations is for an isothermal sphere model 
with Rg = 60 pc. For this class of model there is a large parameter space with Ri, between 10 
and 20 parsecs and X-ray luminosities less than 4.8 x 10^^ erg s~^. For example, for a = 0.2 
and Cf = 0.5, Rb = 16.5 pc, L,j. = 1.2 x 10^^ erg s~^, and P5 = 1.1 x 10"^° dynes cm~^. 
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X-RAY SURFACE BRIGHTNESS PROFILES 



The morphology of the X-ray emission provides a diagnostic to test our models. Match- 
ing up the X-ray emission (IHamaguchi et al.l l2007al ) and the visual images of the Carina 
nebula JSmith fc BrookslboOTbh . Fig. [HI it can be seen that the X-ray surface brightness in- 
creases near the edge of the bubble, the latter seen as a silhouette in the visibl e image. This 



strong ly suggests that a significant fraction of the diffuse X-ray emission seen by lSeward et al. 
( I1979I ) was produced by gas evaporating from the bubble wall. 



This sur f ace br ightness profile is not predicted by either the IChevalier fc Cleggl (119851 ) or 
Castor et al.l (119751 ) models; the former predicts that surface brightness falls monotonically 
with increasing r, while the latter predicts a nearly flat surface brightness profile out to 
r = Rb. For an example of the latter, consider one run with a = 1/3, Cf = 0.5, an 
isothermal sphere ISM with Rq = 60 pc, and stellar wind, radiation pressure, and gravity all 
included. At an age of 3 Myrs we find a bubble radius of 34 pc and a total X-ray luminosity 
of 1.3 X lO^^'erg s~^, roughly consistent with the observed bubble. The cumulative X-ray 
emission as a function of the 3D (not projected) radius is shown in Fig. M^a); the X-ray 
profile with 2D (projected) radius is shown in Fig. [9](b). The calculated profiles show an 
al most flat X-ray profile with a sharp edge at the bub ble boundary. This con trasts with 
the Seward et al. ( 1979 ) and more recent XMM images (jHamaguchi et ah 2007a), which are 
much more complex, with an elevated surface brightness near the bubble walls as noted 
above. 



5.1. Globule evaporation 

To estimate the emission from the interface between the hot stellar wind material and 
the molecular gas around the bubble, we first look at a simplified case of a spherical globule 
of cold gas embedded within the hot gas of the bubble interior, the latter having temper- 
ature Tf and number density rij at large distances from the col d globule. The cold ga s 



will be heated by conduction and evaporate, as discussed in, e.g., ICowie &: McKed ( 119771 ). 
The theory makes three assumptions; firstly the globule is large enough to give a roughly 
time-independent solution for the mass loss rate, secondly that the thermal conduction is 
unsaturated, and thirdly that the fiow has a moderate Mach number {Ai"^ << 5). The con- 
duction is unsaturated if the thermal mean free path for electrons is much shorter than the 
temperature scale height. For the Carina nebula, with ~ 0.1 cm~^ and Tj. = 6 x 10^ K, 
the conductivity is approaching the saturated limit for globule radius i?i ~ 1 pc. Assuming 
we are in the unsaturated case, the mass loss rate is given by the classical equation (same 
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as equation 2) 

. lQ'iTiiK{T)R^ 

where is the mean mass per particle, K{T) is the conductivity at Tf, and Ri is the radius 
of the globule. The dimensionless form of the equation of motion is 

^ln^^_ 2(6-5,) + ^^-l 

where is the mach number of the flow and y = r/Ri. We assume that the flow will shock 
somewhere close to the sonic point and thus still satisfy << 5 at all radii. 

The solution of the energy equation gives the temperature profile 

T{y) = Tj{l-y-'f\ (15) 

assuming the globule temperature T << Tf. The local sound speed can be calculated from 
the temperature profile and thus the velocity deduced from the local Mach number. The 
local density is then 

m 

Any^RfCsM 

The run of temperature and density can be combined with the expression for Bremsstrahlung 
X-ray emission to give the expected X-ray luminosity and surface brightness for an evapo- 
rating globule. The intensity along a line of sight with an impact parameter 6 to a globule 
is proportional to 

dzn\z)T'/\z)e-'''/'''' ^ J dz'^^^^^e'^^/'^^ ~ (^^^ (17) 

where z is measured along the line of sight, 2; = in the plane of the sky at the location of 
the globule, and r = \/W+z^. Once T(r) approaches hu/k, the surface brightness depends 
primarily on b, falling rapidly with increasing b as long as -Ri < b. Doing the integral 
numerically, setting Tf = 6 x 10^ and -Ri = 1 pc we get m = 6 x lO^^g s~^ and the X-ray 
emission profile shown in Fig. [TUl the dotted line shows the b~^ scaling. In words, a spherical 
globule should be surrounded by a fairly sharply defined X-ray halo. 

As a check, the dependence of the total X-ray luminosity and mass loss rate on Tf and 
Ri were investigated. We varied the external temperature for a globule of i?i = 1 pc; in 
a separate set of integrations at fixed Tf = lO^K we varied the radius of the globule Ri. 
These simulations show that M oc T^-^, oc T, and M oc T as expected analytically for an 
unsaturated conductivity. 
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A very crude calculation of the X-ray emission produced by the conductive evaporation 
of a collection of globules forming the bubble wall suggests that most or all of the X-ray 
emission could arise in this manner: 



47rr^A^n^(r)cir. (18) 



Using eqn. f|T6|) . 



yAnCsTripJ Ri \lpcj \3 X 10 "'■^ergs ^cm-*/ 

If the bubble wall (or about half of it) were tiled with /2i ~ 3 pc globules, evaporation from 
the wall could explain the observed X-ray emission. 

In Carina the situation is much more complicated than this toy model. The density of 
the hot gas is manifestly not uniform, and in addition we have seen that hot gas apparently 
escapes rapidly from the bubble interior. We are assuming that this implies that the gas 
density decreases rapidly away from the bubble wall, on a scale comparable to the size of 
the wall fragments, which are typically much smaller than the radius of the bubble. Further 
work along these lines is clearly needed. 

We have extracted sur face brightness profiles al ong rays emanating from t] Carina using 



the XMM observations of (IHamaguchi et al.ll2007al ). The southern part of the XMM image 



is shown in Fig. [TTl which also shows the three rays along which the surface brightness was 
calculated for Fig. [T21 These profiles were produced by interpolating through the pixels 
using IDL and then smoothing. All three lines were interpolated through 2554 points and 
then smoothed by a box car average with a width of 60 points. 

Emission was converted from counts to flux using the PIMMS package. We assumed a 



spectrum of two black bodies as found by IHamaguchi et al.l (l2007al ). 



The 700-1300eV band shows the greatest contrast, from little emission in the bubble 
interior (but away from rj Carina) to substantial emission near the optically dark lanes. The 
400-700eV band also shows some edge effects, especially along lines Two and Three. The 
high energy band, 2000 — 7000 eV, shows essentially no enhanced emission near the bubble 
walls. 

In the 700-1300eV band line One shows a decrease in surface brightness from rj Carina 
out to 200 arcsec, a fairly flat proflle between 200 and 650 arcsec, followed by a decrease out 
to the edge of the image. Line Two shows a similar flat proflle between 150 and 400 arcsec, 
decreases to a minimum surface brightness at 500 arcsec, then increases to a peak at about 
900 arcsec. The proflle along line Three is similar to line two but with a less pronounced 
minimum, and a peak closer to 1] Carina (at 750 arcsec). The peak along ray Three is wider 
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than that along ray Two as it makes a larger angle with the bubble edge normal. These 
profiles quantify the visual impression that the bulk of the bubble interior has a fairly flat 
X-ray surface brightness but that there is strong emission from the edges. The edge profiles 
are of significantly different shape to that expected from a Castor bubble, Fig. M^a) and are 
more curved than would be expected from globule evaporation, Fig. [TOl 



6. DISCUSSION 



We have noted three discrepancies between the Castor et al. scenario and observed 
bubbles: the predicted bubble radii are too large at a given cluster or stellar age, the predicted 
X-ray luminosities are too high by a factor of 100 or more, and concomitantly, the predicted 
pressure of the hot gas is too large by a factor of 10 or more. Finally, we have noted that the 
radio free-free emission is consistent with a filling factor of H II gas that is of order unity, 
suggesting that the pressure in the bubbles is controlled by the H II gas rather than that of 
the shocked stellar wind gas. 

A number of aut hors have noted the small observed - R/,, roughly by a factor of 5, com- 



pared to predictions (IDorland. Montmerle. & Doom 


1986; 


Dorland & Montmerle 


1987; 


Oey 


1996; 


Rauw et al. 


2002; 


Dunne et al. 


2003; 


Smith. Struck. & Nowak 


2005 


). These authors 



have also noted the deficit of X-ray luminosity or hot gas pressure compared to predictions. 
Suggested resolutions involve one or more of the following: lower stellar wind luminosities, 
mass loss from the bubble, or energy loss from the bubble, e.g., from super-Spitzer effective 
conductivities, or highly efficient mass loading so that the shocked stellar wind cools by 
conduction below X-ray temperatures. 



Castor et al.l (Il975l ) and IWeaver et al.l (119771 ) make three assumptions; first, that the 



energy deposited by the stellar winds (or 5/11 of it) is stored in the hot gas inside the bubble, 
second, that the hot gas mass is controlled by conduction (at the Spitzer rate), and third, 
that the surrounding ISM is uniform (represented by their constant density no). Applying 
the model requires knowledge of L^, no, and the age of the system. The predicted bubble 
radii depend only weakly on and uq, while the ages of star clusters are reasonably well 
constrained (to within a factor of 2 — 4). Since the bubble radii are universally overestimated, 
at least one assumption must be incorrect. 

We can acknowledge the fact that the dynamics of the model are wrong, but still test the 
first assumption, that the energy deposited by the stellar winds is stored inside the bubble. 
We will also assume that conduction drives gas from the bubble wall into the interior until 
Tb ~ 10^ — 10^ K; this assumption is bolstered by the observed gas temperatures, which 
cluster around 6 x 10^ K. 
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The argument is that the time integrated wind luminosity deposits an energy E = 
(5/ll)L^r in the bubble interior, and that the corresponding pressure is 2E/3V, where 
the volume is assumed to be a sphere with the observed bubble radius. Together with the 
assumption that the temperature is T ^ 5 x 10^ K, we can calculate the mean density and 
hence the X-ray luminosity. 

Both X-ray and other observations give estimates of the bubble radius. Star counts 
(to a given limiting magnitude) and radio observations give estimates of L, the bolometric 
stellar luminosity, or Q. The stellar wind luminosity can then be estimated from either L or 
Q, with the caveats mentioned above. Finally, from the presence of Wolf-Rayet stars and/or 
SN remnants, we have estimates of the cluster age. In terms of these quantities, the X-ray 
luminosity is 



3 X 10^^^ 



4 X 10^^ erg s^^ s^^ 



20 pc 



6 X lO'^KA^ /3.6 X lO^yr^ ^ 



T 



T 



(20) 



where we have assumed an X-ray cooling rate A^; ^ 3 x 10^^'^^ erg s~^ cm^, e.g., IChu et al. 
( 1l995l ). ^ is the metallicity relative to solar, and r is the age of the cluster. We have scaled 
to values appropriate to Carina. As previously noted, this is 10,000 times too high. 

In addition to Carina, we have calculated this predicted X-ray lumino sity for the LMC 
H II r egions detected in X-rays by the Einstein satellite, and discussed in IChu fc Mac Low 
( ll990l ). We estimate the bolometric luminos ity of the H II regions in two ways. First, we use 
the star counts from lLucke fc Hodgd (1l970l ): these authors report star counts for stars with 
niv < 14.7, about My = —3.85 at the distance of the LMC. Second, we use the observed 



radio free-free flux from iMcGee. Brooks, fc Batchelorl (119721 ) to estimate Q and thence L; 
the results from the two methods are consistent, within a factor of 3 or so. As in the case 
of Milky Way bubbles, this indicates that the filling factor of the H II gas is of order unity, 
and hence that the pressure in these bubbles is not set by the X-ray emitting gas. From the 
bolometric luminosity we then estimate L^. Using these values, eqn. ( l20l) over-predicts the 
X-ray luminosity, by factors ranging from a few to five hundred. 

Thus the picture for the LMC is consistent with our finding for Carina. Confusingly, 
other authors have found an excess of observed X-ray emiss i on compared to predictions o f 
the Castor et al. model Jchu fc Mac Lowlll99d fchu et allEggsl : IWang fc Helfandl E99l[ ) . 
Why should this be? 



Chu fc Mac Lowl fll990l ): IChu et all (Il995f ) and lWang fc Helfandl (Il991f l all use the seal 



ing laws of Castor et al. to predict the X-ray luminosity of bubbles in the LMC. The latter 
two give expression for that appear to contain only directly observed quantities, such as 
the bubble radius and the bubble expansion velocity. However, the coefficient relating to 
the observed quantities depends on the theoretical estimate for the bubble radius and central 
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bubble density. This dependence results in an estimate for the stellar wind luminosity that 
is more than a factor of 10 smaller than the wind luminosity as estimated from either simple 
star counts, or from the free-free radio emission. 



In the case of IChu fc Mac Lowl (Il990l ). this can be seen explicitly in their discussion of 



N51D, where they infer L 
star counts of 



Lucke fc Hodgd ( 



from iMcGee. Brooks, fc Batchelorl (Il972l ) we find 
value, and the observed R}, = 48 pc eqn. fl20|) yields L 
5 larger than the observed Lx = 3.3 x 10^^ erg s~^. 



..08 X 10^^ erg s ^ in their appendix B section IV. From the 
1970) we estimate 

1 



2 X 10"^" erg s ^; from the radio flux 



^38 

6 X 10^^ 



erg s 
1.6 X 10^6 



erg s 



Using the lower 
^, about a factor 



Since their estimates of are smaller than the observed X-ray luminosities, IChu fc Mac Low 



(1l990l ) suggest that the X-ray emission of the LMC objects is provided by supernovae that 
have gone off near the bubble wall. We argue that, using the estimate given here, recourse to 
supernovae is unnecessary. In addition, the morphology of the X- ray emission in Carina d oes 
not support this picture, nor does that of the objects studied by lChu fc Mac Lowl (Il990l ) in 
the LMC. In Carina, as stated above, the emission is slightly enhanced near the bubble walls, 
but not at a particular point, as would be expected following a supernova. 

Similarly, the X-ray images presented by iTownsley et al.l (120061) of N 157 (30 Do- 
radus), especially their Figures 14 and 15, generally show rather diffuse emission in the 
cavities surrounded by cold dusty gas, reminiscent of what is seen in Carina. However, as 



Townsley et al.l (120061 ) also remark, in some cases the emission comes from the interior of 



the cavity, while in others it comes from edge brightened regions with distinct central voids, 
similar to Carina. Tw o exceptions may be the regions referred to as numbers 5 and 9, which 



Townsley et al.l (120061) suggest rn ay be the result of supernovae interacting with the bubble 



wall, a la lChu fc Mac Lowl (119901 1. 



6.1. Dynamics Controlled by H II Gas 

As we noted in section ^ the H II gas pressure in Carina, as estimated from the radio 
free-free emission, is similar to the pressure in the X-ray emitting gas; both are similar to 
the pressure w here the outflows frona photo -dissociation regions shock. This was pointed 



out earlier by iDorland &: Montmerld (119871 ). Since the sound speed of the cooler gas is 
comparable to the bubble expansion velocities observed, this is not surprising. 

What is more surprising is that the free-free luminosity is consistent with a filling factor 
for the H II gas that is of order unity (see eqn. 12|) a fact apparently not noted before. This 
implies that the gas pressure is set by the H II gas rather than the hot gas. It suggests that 
the shocked stellar winds escape from the partially confining bubble at the sound speed of 
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the hot gas, as long as the pressure in the hot gas exceeds the H II gas pressure in the region. 
When enough hot gas has escaped that ^ Prii, the escape will slow, as the hot gas is 
impeded by the cooler gas. This sets up a rough pressure balance between the hot and H II 
gas. 



Borland. Montmerle. &: Doom! (1l986l ) and iDorland fc Montmerld (119871 ) consider the 



loss of wind energy (by an enhanced conduction by a non-Maxwe llian tail of high energy 
electrons, according to the latter). iDorland fc Montmerld (119871 ) note that the thermal 
energy in hot gas in the Rosette as well as the Carina nebula cavities is well below that 
input by stellar winds over the age of the respective star clusters. These authors suggest 
that the excess energy is conducted away from the hot gas to surrounding cold gas, rather 
than by simply flowing out through holes in the bubble wall as suggested here and elsewhere. 
This would alleviate the dynamical problem (the predicted bubble radius would be smaller 
than in the Castor et al. theory) and the X-ray luminosity problem. 

In a similar vein, iMcKee. Van Buren. fc Lazarefj (119841 ) consider the evolution of wind 
blown bubbles in conjunction with H II regions. Like us, they conclude that the pressure 
and hence the dynamics of the bubble are controlled not by the 10^ K gas, but by the H II 
gas. They suggest that the energy density of the shocked wind is dissipated in situ, rather 
than leaving the bubble. In the case of a high wind luminosity (as in Carina), they suggest 
that the hot gas will engulf a large number of globules, producing a large mass loading and 
resultant catastrophic radiative cooling. 

Examination of the Hubble Carina mosaic (jSmith &: Brooksll2007bl ). which has a spatial 
resolution of ~ 10~^ pc, does not reveal a large population of evaporating globules in the 
interior of the bubble. If the bubble wall consists of a large number of sub-resolution globules 
it might be possible to dissipate the wind luminosity; the bulk of the wind luminosity would 
be radiated near or just beyond the Lyman edge, or possibly as optical light. Radiation 
beyond the Lyman edge would be hard to detect due to absorption along the line of sight 
to Earth, while optical emission would be hard to distinguish from emission associated 
with ionizing radiation from the star clusters. However, preliminary calculations of the 
emission region associated with an evaporating globule suggest that a few percent of the 
luminosity would emerge in the X-ray band; only if the globules are ~ lOOAU in radius 
is the fraction of emission emerging in the X-ray band small enough to escape detection. 
The conduction around such globules is highly saturated; in fact, the mean free path for 
electrons to thermalize is A ~ 0.1 pc. Hence, the globules would have to be separated by 
many times their own radii to avoid an X-ray emission fraction larger than observed. These 
small globules are surrounded by a saturated conduction region roughly ten times their own 
radius; if these saturated conduction regions overlap they will radiate X-rays like a single, 
much larger globule. We defer further discussion of this possibility to a later publication. 
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Neither lDorland. Montmerle. fc Dooml (119861 ) nor lMcKee. Van Buren. &: Lazarefa (Il984j ) 
explain why Phii ~ ^/^QJA^k^/x^ (where the latter assumes a filling factor near unity). It is 
not clear why the energy-loss rate of the hot gas would tune itself so that the filling factor of 
both warm and hot gas is near unity. A large cooling luminosity would lead to an essentially 
complete absence of hot gas, while a low cooling luminosity would result in too high a hot 
gas pressure. Unlike a simple mass loss scenario, an energy loss scenario requires some fine 
tuning to match the observations. 

Finally, we note that our leaky bubble models, which predict smaller bubble radii and 
gas pressures, have implicat ions for the predictio n of spectral energy distributions of starburst 
galaxies. As pointed out by lDopita et al.l (120051 ). using a standard Castor et al. bubble model 
results in ionization parameters too small (by a factor of 10 or more) than that inferred from 
spectral energy distributions of such galaxies. Both the smaller radii and the lower gas 
pressure (hence gas density) will increase the ionization parameter associated with bubbles 
surrounding massive star clusters. 



6.2. Ballistic Bubbles 

In the Castor et al. scenario, if the shocked stellar wind gas from a star or star cluster 
was confined to a bubble of volume V , the pressure of the hot gas would be 

where the energy in the bubble Ei, = {5/11)Lj^t and r is the age of the star cluster. As 
already mentioned, this pressure is orders of magnitude larger than the observed pressure in 
Carina, so the hot gas is not currently exerting the predicted pressure on the surrounding 
dense gas. 

However, it is possible that in the past the hot gas was confined. The force exerted by 
the gas on the bubble would have been much larger than that exerted by the H II gas or by 
radiation, resulting in a nearly impulsive expansion of the bubble. After this early impulsive 
phase, the hot gas could then leak out of the bubble, leaving the bubble to expand in a nearly 
ballistic manner. This could happen, for example, if the density of the trapped hot gas grew 
large enough tha t the radiative cooling time scale became shorter than the expansion time 



scale ( 1Shulllll980l ). We argue that this is not the case. 



First, high re s olutio n radio images of ultracompact H II (UCHII) regions {R ~ 0.1 pc) 



by iDe Free et al.l (120051 ) show that most are either shell-like (see Fig. [T3l) or cometary (a 
significant fraction are bipolar). The former two would allow hot gas to leak out of the 
UCHII region. Observations show that the expansion velocities are typically ~ 10 km s^^. 
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e.g., (jChakraborty fc Anandarad 119971 . Il999l ). The expanding shell will sweep up mass; if 
the motion is ballistic the velocity will decrease and the bubble will never reach sizes of order 
lOpc. 

Second, there are U CHII regions that emit diffuse X-rays, e.g., in the Welch ring of W49 



(|Tsuiimoto et al.ll2006l ) and in the Hourglass region of the Lagoon nebula, M8 (IRauw et al. 



2OO2I ) . In both cases the X-ray luminosity is well below that expected under the assumption 
that the shocked wind gas is confined in a bubble; it is also well below L^. We discuss each 



source in turn. 



The W49 ultrac ompact {r ^ 0.15 pc, D = 11.4 kpc) source Gx has ^ 3 x lO'^^erg 
s~^ and T 7keV (iTsujimoto et al.ll2006l ). much hotter than the plasma in Carina. The 
high temperature of the X-ray gas is difficult to reconcile with the Castor et al. model; 
such hot gas is more likely to occur in a wind-wind collision. If this is the case, then the 
luminosity of any truly diffuse, shocked, partially confined gas is even smaller than the 
value given here. The column toward the source is S ^ 0.8gcm~^, suggesting a dynamical 
pressure P = t^GT? 10~^ dynes cm~^, and a mean cold gas density p ^ 2 x lO^^'^gcm"^. 
This is consistent with their estimated hot gas pressure based on the emission measure 
{ux ~ 10cm~^, Px ~ 10~^ dynes cm~^). They suggest the driving source is one or two 05 or 
06 stars, with a luminosity in the range 2 — 6 x IO^Lq, consiste nt with the ion i zing p hoton 
luminosity Q ~ 10^^ s~^ for each of the five G sources listed by iDe Free et al.l (119971 ). The 
corresponding wind luminosity is ~ 500Lq. 

The age of the source can be estimated in two ways. First, using the mean density, the 
dynamical time is r^yn = l/^/Gp 3 x 10^^ s. Second, using the force exerted by the hot. 
X-ray emitting gas; as just noted, the latter has a much lower density than the absorbing 
gas, but a similar pressure. The outward acceleration is P^/S ~ 10^^ cms"^. The dynamical 
time is then r^y„ 2 x 10^^ s, roughly the same as the estimate based on the self-gravity of 
the cold gas. 

Assuming that the stellar wind is trapped in a Castor style bubble, the energy in the 
shocked gas is E = ■^[LujT ^ 4 x 10^^ erg s~^, and the pressure P = 2 x 10~^ dynes cm~^, 
two orders of magnitude higher than observed (and with an associated X-ray luminosity four 
orders of magnitude too high). Since this pressure is so much higher than the dynamical 
pressure, and the sound crossing time of the hot gas so much smaller (by a factor of about 
1000), the simplest interpretation is simply that the shoc ked stellar win d gas e scapes, as in 
Carina. The radio images of the G sources in Figure 2 of iDe Free et al.l (120001 ) show a very 
incomplete ring, strengthening the argument for rapid escape. 

The second case is the ultracompact H II region G 5.97 -1.17 associated with the Hour- 
glass nebulae in M 8. It is believed to b e powered by the 07 V star Herschel 36; using 
values for M^, from iRepolust et all (120041 ) (not allowing for a clumped wind), we estimate 
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Lyj = 2 X 10^^ erg s~^, while iRauw et al.l (120021 ) use a wind luminosity a factor of 3 smaller. 
They estimate a dynamical age of 2 x lO^yr based on the size of the Hourglass (they use 
R = 0.2 pc ) and the observed expansion velocity of 10 km s~^ ( Chakraborty fc Anandarao 
1997l . ll999[ ). Using their values for L2, t^™; and Rb, we find a predicted hot gas pressure 



of P = 4 X 10"''' and Ux ~ 200 cm~^. The predicted = 7 x 10^^ erg s~^, larg e r than that 
observed by a factor of 1000. Using the full Castor et al. model, iRauw et al.l (120021 ) note 
that the predicted bubble radius is larger than that observed (but only by a factor of 2) and 
that the predicted X-ray flux is ~ 10^^ erg s~^. This smaller than our prediction by a 
factor of 7, consistent with their use of the (two times too large) Castor et al. radius. Their 
predicted X-ray flux is still larger than observed by a factor of 175. 

Radio observations by [Turner et al.l (119741 ) find a radius r ~ 0.1 pc (half that sug- 
gested by Rauw et al.) and a flux at 11cm of 0.54 Jansky. The H II gas pressure is 
1.5 X 10^^ dynes cm~^ from the observed radio flux, and 1.5 x 10^^ dynes cm~^ using Q 
for an O 7 V star. This should be compared to the pressure of the X-ray emitting gas 
Pj, = 8 X 10~^ dynes cm~^, and to the radiation pressure L/Anr'^c = 1.6 x 10~® dynes cm~^. 
All these pressures are at least a factor 10 below the pressure found under the assumption 
that the wind luminosity is trapped in the bubble. We conclude that, as for Gx, the most 
likely explanation of the low observed X-ray flux and the small bubble radius is that the 
shocked stellar gas simply escapes from the bubble interior. 



7. CONCLUSIONS 

We use both the Castor et al. theory and the Chevalier and Clegg theory for bubble 
evolution around massive star clusters to try to understand observations of Milky Way and 
LMC bubbles. The Castor theory over-predicts X-ray luminosity in the Carina bubble by 
a factor of 60 and expands too fast by a factor of 4. In contrast, the Chevalier and Clegg 
model under predicts the X-ray luminosity by a factor of 20; since it is a steady state it has 
no expansion time. 

These results suggest there is a partially confined hot, homogeneous interior but there 
are also holes in the shell through which the hot gas escapes like a free flowing wind. We 
constructed a model in which the confining shell covered only a fraction Cf of the sky as seen 
from the cluster, and allowed hot wind gas to escape at its sound speed from the uncovered 
portion of the shell. As an example one model consistent with observations of Carina is 
an isothermal sphere CMC model with Rq = 60 pc, a = 0.2, and Cf = 0.5, Rb = 16.5 pc, 
Lj, = 1.2 X 10^^ erg s~^, and Pb = 1.1 x lO^^'' dynes cm~^. This compromise is the best 
way to fit the many observations of the Carina nebula and is likely the best solution for all 
bubbles. Our models suggest that the exact size, distribution and evolution of holes does 
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not matter much, it is the weighted average of the fraction of total area that matters. 

The bottom hne is that the dynamics of X-ray bubbles are not strongly affected by 
the presence of luminous stellar winds; rather, they are controlled by a combination of the 
self- gravity of the surrounding G MC gas, radiat i on pr essur e, and H II gas pressure. This 
conc lusion is in contrast to that of ICastor et al. ( 1975 ) and ShuU ( 1980 ). but in agreement 
with iMcKee. Van Buren. fc LazarefB (119841 ) . We argue that the X-ray fluxes of the bubbles 
are always below those predicted by the Castor et al. picture, properly interpreted, i.e., using 
the ob s erved rather than predi c ted b ubble radii. As noted by lDorland. Montmerle. fc Doo 
(119861 ): iDorland fc Montmerld (119871 ) , the X-ray observa tions sh ow that the pressure in the 
hot ga s is well below that predicted bylCastor et al.l (|l975l): unlike iDorland. Montmerle. fc Doom 
(119861 ) ; IDorland fc Montmerld (119871 ) and iMcKee. Van Buren. fc Lazarefj (119841 ) , we argue 
that this is a result of leakage of hot gas from the bubble, rather than leakage of energy in 
the form of radiation. 

The results of the models in this paper are dependent upon the ISM model used, both the 
run of mean density, and the fluctuations in the surface density in different directions from 
the bubble center. According to both simulations and observations, the latter is distributed 
log-normally, meaning that in many directions the surface density will be far below the 
mean surface density. We suggest that hot gas escapes along these directions, which in our 
simplified models we treat as holes in the bubble wall. 

Full hydrodynamical models are needed to check that the semi-analytical solutions 
shown in this paper are correct and to understand both the early evolution of bubbles and 
the development of holes in the shell. 
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A. CHEVALIER & CLEGG BUBBLES 

Chevalier & Clegg assume a spherically symmetric wind with mass input and energy 
input E = Luj. For a smooth transition from subsonic flow at the centre to supersonic flow 
at large distances the Mach number must equal one at some radius, which Chevalier & Clegg 
take to be r = Rd. 

The flow is described by 

i|:(p«r^)=0, (Al) 
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where u is the wind velocity, r is the radial co-ordinate, p is the density, P is the pressure, 
and 7, the adiabatic index, = 5/3 throught this paper. For r < Rd, (f) = M^/V, $ = E/V, 
V = 47ri?y3; for r > Ra, = $ = 0. Equations have the solution 



for r < Rci, and 
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for r > Rci, where is the Mach number. 

The flow can be integrated in two sections, r < R^ and r > R^ using the outer boundary 
conditions of the first as the inner boundary of the second. Thus, it is possible to follow 
the evolution of the wind and numerically integrate the expected X-ray luminosity from the 
cluster edge out to arbitrarily large distances.. 

Initially a small radius (r << Rd) and small initial mach number, A4o are assumed. 
The true Mach number is then found iteratively to 10 decimal places for that radius 
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The calculated Mach number was then combined with the observed mass loss rate and stellar 



wind luminosity for Tr 16 (Table 1) to find the initial conditions (from I Chevalier fc Clegg 
f |l985h ) 
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From these initial conditions equations flAll) - flA3l) can be integrated to the cluster radius. 
The integration continues with the evolution equations (without the source terms outside 
the cluster) and the appropriate Mach number iteration, 

2 1_ 



RciJ \l + 3/M, 



B. CASTOR et al. BUBBLES 



For a bubble in the snow plough stage, Rp is large enough tha t the flow from c to b can 
be considered stationary and plane parallel. IWeaver et al.l (119771 ) considered the effects of 
relaxing this assumption and found the results are only changed by <10%, much less than 
the observational errors fo r, e.g., the bubble radius in Carina. For a stationary plane parallel 
flow at constant pressure (jPenston fc Brownlll970l ) 



dH _ d 
^ dz dz 



K(T)f 

dz 



pHt), 



(Bl) 



where T is temperature, z = Rb — r, pv = Mh/AnRf, is the mass in region b, if = bk^T /2p 
the specific ent halpy, K{T) = CT^^"^ is the thermal conductivity, and C = 1.2 x 10~^ erg 
cm~^ s^^ K''/^ ( ISpitzeii Il962l ). Integrating equation (]Bip from the b-c boundary: z = 0, 
T ~ to the centre: z = Rb, T = Tb gives the mass evaporating from c to b as 



Mb 



'25 h 



K{T)Rb. 



:b2) 



To get this solution we have assumed the cooling term, pA{T) is negligible in the transition 
zone between b and c. 

Energy is constantly added into region b by the stellar winds with the total luminosity 
Lw and the dominant energy loss is through the expansion of the bubble. The energy in 
region b obeys 

Eb = L^- 4:TTRlPbRb, (B3) 
where Pb is the pressure in region b and 

4 



nRtPb 



-Eb 



(B4) 



When considering radiative losses from region b it is necessary to follow the evolution 
of region a (the hypersonic stellar wind region) as eventually region b will cool enough to 
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collapse. Therefore, the equations need to be altered to account for the fact that the size of 
region a is not negligible. In that case, the energy in region b is given by 



E, = 2-kP,{RI - Rl) 



(B5) 



where Ra is the radius of region a. Ra can be calculated from pressure balance at Ra- The 
pressure from region b on a shell at Ra is 



Eh 



2ttRI-RI- 



(B6) 



The stellar wind luminosity is L^, = so the rate of momentum deposition by the 

stellar wind is 



This then gives a pressure at radius Ra of 
Setting these pressures equal gives 

Ra 
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Radiative losses with i n regi on b can be calculated by reference to the cooling functions in 
Sutherland fc Dopital jlQQ.sh 



L,^nl^7r{Rl~Rl)A{T). 



(BIO) 



Including radiative losses into the Castor et al. evolution equations results in the fol- 
lowing closed set of evolution equations 
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where Vc is the momentum of the shell and Ci = 4.13 x 10 cgs. 
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The initial conditions of the bubble are taken from lWeaver et al.l (119771 ) at time to = 10^ 



years: 
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These give initial conditions close enough to the true values that the numerical solution 
settles down by the cluster boundary and is a valid solution for Rb > Rd- 

Initially the parameters were set up to match those used in IWeaver et al.l (119771 ) and 
results compared to Figure 6 of that paper as a check of our code. Once the code was tested 
parameters for Tr 16 could be used (from Table 1). 
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Fig. 1. — Cartoons of the Castor et al. bubble (i.) and the modified Castor et al. bubble (ii.). 
The Castor et al. bubbles have a four zone structure. Zone a is hypersonic stellar wind; zone b is 
shocked stellar wind and evaporated material from the shell; zone c is the bubble shell; and zone d 
is the ambient ISM. ii. - The structure of the modified Castor et al. bubble, showing leaking from 
the shell. Note the larger size of region a, smaller size of region b and material escaping into region 
d. 
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Fig. 2. — Plots of the basic ICastor et alJ (jl975l ) model evolution against time (a) and Rb (b) for 
Tr 16. These plots show how pressure, number density and total X-ray luminosity change through 
the evolution of the bubble. Pressure and X-ray luminosity have been scaled to fit on the graphs. 
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Fig. 3. — Plot of the steady state IChevalier k. Clegel (|l985l ) model for Tr 16. The number density, 
pressure and X-ray luminosity against radius are plotted. Pressure and X-ray luminosity have been 
scaled to fit on the graphs. Note that we plot the cumulative X-ray luminosity as a function of 
radius. 
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Fig. 4. — Plot of the scaled X-ray luminosity and pressure calculated as a function of radius for 
the two classes of (unmodified) models for Tr 16, as discussed in §4.1[ The Castor model is plotted 
as a function of bubble size whereas the Chevalier model is plotted as a cumulative function of 
radius (due to its steady state nature). The upper and P lines are for the Castor model and the 
lower lines for the Chevalier model. The two dots show the observed pressure (lower) and X-ray 
luminosity (upper). In the Castor et al. model the pressure is roughly constant for r < R^, but this 
constant pressure is a function of the bubble radius, as shown. In the Chevalier &: Clegg model 
there is no bubble, but pressure decreases with increasing radius. 
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Fig. 5. — Graphs for three modified Castor simulations with Cj = 1 (sohd), Cj = 0.65 
(dashed), and Cf = 0.3 (dotted), (a) radius versus time, (b) pressure versus radius, (c) 
X-ray luminosity versus radius, and (d) number density inside the bubble versus radius. 
Observed results are shown as a black dot. See the discussion in §4.21 
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Fig. 6. — Graphs showing the details of modified Castor simulations with different ISM 
models and different shell cover factors Cf. All runs have Rd = 3pc, Rg = 73.86 pc and 
a = 1/3. (a) is Rb versus time, (b) pressure versus Rb, (c) X-ray luminosity versus Rb, and 
(d) rib versus Rb- The solid line is for an isothermal GMC, with Cf = 1, the dotted line is 
an isothermal GMC with Cf = 0.5, the dashed line has constant ISM density with Cf = 1, 
and the dash-dotted line has constant ISM density with Cf = 0.5. Note that (b) and (c) are 
plotted against radius so although some models do fit the observed quantities at that radius 
they do so after too short a time. 
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Fig. 7. — Parameter study results for varying a and Cj for either an isothermal sphere 
ISM and for Rq = 24 pc (a) or Rq = 60 pc (c), or for a constant density ISM for Rq = 24 
(b) or Rq = 60 pc (d) modified Castor models at 3Myrs. The contours show the X-ray 
luminosity (dotted) and the radius (sohd). Gray dashed hnes show the boundary of the 
collapse region. X-ray luminosity is in logarithmic contours in intervals of xlO erg s~^ for 
the homogeneous models and x 10°"^ erg for the isothermal models. Radius contours are 
hnear in intervals of 10 parsecs. The observed properties of Carina are 13 < i?6 < 20 pc and 
La; = 4.8 X 10^^ erg s~^. For the isothermal ISM model there are regions of the parameter 
space where the shell collapses inwards in the simulations. This is due to the swept up mass 
being so large the gravitational force is larger than the sum of the gas pressure force and 
radiation pressure force. 
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Fig. 8. — Image of the Carina nebula in visual (jSmith &: Brooksl l2007al ) and X-ray contours 
(|Hamaguclii et al.ll2007al ). The visual image is a combination of H II, O IV and S II and the X-ray 
image is between 0.7 keV and 1.3 keV. The green outline shows the edge of the X-ray image. See 
the electronic edition of the Journal for a color version of this Figure. 
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Fig. 9. — Graphs of a modified Castor model showing (a) the cumulative X-ray luminosity 
(solid) and X-ray profile (dashed) with 3D radius; (b) shows the X-ray profile with projected 
radius (z); we expect an almost flat X-ray profile with a sharp cutoff at the bubble edge. 
This model is for a = 1/3, Cf = 0.5, isothermal sphere ISM model with Rq = 60 pc. 
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Fig. 10. — The X-ray emission profile against r/Ri for Tf = 10^ and i?i = 1 pc for an evaporating 
globule (solid line) showing a fairly sharply defined X-ray halo. The dotted line shows a gradient 
of -3. 
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Fig. 11. — XMM image of the southern part of the Carina Nebula between 700 eV and 1300 eV 
with contours, rj Carina can be seen along the top (with all three profile lines pointing at it). 
Thick lines are the profile lines used in Fig. [T2j 
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Fig. 12. — Surface brightness along the lines in Fig. [TT] for the three different XMM bands. The 
solid line is for profile one, the dashed line for profile two and dotted line for profile three. Surface 



bright ness is given in erg cm ^s i, with an assumed Nh = 1-8 x lO^i cm ^ ( Hamaguchi et al 



2007al ) : projected distances from t] Carina are given in arcseconds. Recall that at a distance of 



2.3 kpc, 1000 arcseconds is about 11 pc. 
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Fig. 13. — VLA 7mm observation of W49A/D from lDe Pree et al.l (|2005l ). Beam size shown lower 
left is about 500 AU. The shell has a radius of ~ 0.03 pc and the column density through the shell 
is significantly variable. 
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Number of stars > B2 


72 


Total bolometric luminosity 


1.7 X 10^ Lq 


Total ionizing flux (Q) 


5.9 X 10^° s-i 


Total mass loss rate Myj 


1.1 X 10-3 Mo/year 


Total stellar wind luminosity 


6.7 X 10^ Lq 



Table 1: Properties of Tr 16 in the Carina Nebula using all stars of spectral type B2 or earlier 
assuming standard mass loss rates (jRepolust et alll2004l ). From ISmithI (|2006al ). 



T/K 


Ld 1 erg s"^ 


Qe 


Md/M^ 


35 


2.93 X 10^° 


8.72 X 10-^ 


6 610 


80 


1.19 X 10^° 


4.71 X 10-4 


18.2 


220 


4.24 X 10^9 


3.72 X 10-3 


0.014 



Table 2: Dust masses of the three components seen in ISmith &: Brooks! (j2007al ) with better dust 
modeling, compared to Smith &: Brooksl ( 2007a ) total dust mass of 9.6 X lO^M© 



